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Quantum fluctuations stabilize Skyrmion textures 
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Here we show that the zero point energy associated to the quantum spin fluctuations of a non- 
collinear spin texture produce Casimir-like magnetic fields. We study the effect of these Casimir 
fields on the topologically protected non-collinear spin textures known as skyrmions. We calculate 
the zero point energy, to the lowest order in the spin wave expansion, in a Heisenberg model 
with Dzyalonshinkii-Moriya interactions chosen so that the classical ground state displays skyrmion 
textures, that disappear upon application of a strong enough magnetic field. Our calculations show 
that the Casimir magnetic field contributes a 10 percent of the total Zeeman energy necessary to 
delete the skyrmion texture with an applied field, and is thereby an observable effect. 
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The properties of many magnetic materials with bro¬ 
ken spin symmetries are correctly described in terms of 
magnetic moments that have a well defined orientation 
and are, in this sense, classical. Yet quantum fluctua¬ 
tions, that ultimately arise from the fact that the differ¬ 
ent projections of the spin operator do not commute with 
each other, are invariably present. The standard[TJ [2] 
description of spin excitations of these broken symmetry 
magnetic materials in terms of quantum spin waves im¬ 
plies that, even in the ground state, there is a zero-point 
(ZP) energy ultimately associated to fluctuations of the 
magnetization field [1] |2 . This ZP energy vanishes in the 
case of collinear ferromagnets, but is non-zero in general. 

In the case of quantum electrodynamics, the vacuum 
energy, Ezp , associated to the photon zero point motion 
stored between two parallel plates depends on the dis¬ 
tance between L between them, resulting on an effective 
Casimir force[3] F = — 6S S ^ P . In this letter we show that, 
analogously, the vacuum energy associated to the quan¬ 
tum spin fluctuations of a given spin texture depends 
on an ensemble of classical variables encoded in the lo¬ 
cal spin orientation resulting in an effective Casimir 
field Ti = — 6S Z r • Here we show that proper account of 

O u 2 't 

those fluctuations leads to a reduction in the energy of 
skyrmion textures rendering them stable in regions where 
a purely classical analysis predicts them to be unstable. 

Magnetic Skyrmions are topologically protected spin 
structures [4]. Observed recently, both in chiral 
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magnets [5MTT] and in engineered surfaces [17, 18], they 
have received attention for potential applications in spin- 
tronics because it is possible to control their position with 
very low current densities [12;. In addition skyrmions dis¬ 
play several features that are convenient from the view¬ 
point of potential applications, such as dimensions within 
the nanometric scale, topological protection, and high 
mobility. They have been observed in bulk magnets 

MnSi[5H7], Ftei_ a Co x Si0 U US E], Mm^Fe.GeHS] 

and FeGe[T6] using neutron scattering and Lorentz trans¬ 
mission electron microscopy. As it was reported, in these 
systems an external magnetic field can induce skyrmions 
with diameters of about a few tens of nanometers. The 
inclusion of spin transfer torques, as it is shown by nu¬ 
merical simulations, can be used to nucleate and manip¬ 
ulate isolated skyrmions fl9l 120] . In addition, two dimen¬ 
sional atomic-scale magnetic nanoskyrmion lattices have 
also been observed by means of spin polarized scanning 
tunneling microscopy in a monoatomic layer of Fe atoms 
on top of Ir(lll) surface [17] • In this case, skyrmions 
arrays with a length scale of only one nanometer are sta¬ 
bilized by the interplay between Dzyaloshinskii-Moriya 
interaction and the breaking of inversion symmetry at 
surfaces and interfaces. 

Previous work [21] has addressed the study of 
skyrmions spin excitations in a semiclassical approxima¬ 
tion, for rather large systems with up to 10 5 localized 
moments. An effective theory of magnon-skyrmion scat¬ 
tering in a continuum approach has been recently car¬ 
ried out by Schiitte et al [22]. Here we treat smaller 
skyrmions, using a complementary approach that per¬ 
mits to study quantum fluctuations associated to the 
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excitations relevant for atomic scale skyrmions [17] for 
which the effects of quantum fluctuations of the magne¬ 
tization cannot be neglected at the outset. 

The non-collinear spin alignment between first neigh¬ 
bors is promoted by the competition between the stan¬ 
dard Heisenberg interaction ( — JSi • Sj) and the an¬ 
tisymmetric Dzyaloshinskii-Moriya interaction (DMI) 
term Dij • Si x Sj, which is only present in non- 
centrosymmetric systems. In particular, this competi¬ 
tion leads to helical spin textures with fixed chirality. 
Application of a magnetic field h along the z direction 
(perpendicular to the system surface and then to Dij for 
any pair i,j ) breaks the up-down symmetry in the sys¬ 
tem resulting, for a range of fields, in an abrupt change in 
the ground state. In this manner, the external magnetic 
field and a uniaxial anisotropy can stabilize an isolated 
Skyrmion on a ferromagnetic system with DMI. A mini¬ 
mal model that captures the essence of the physics just 
described is given by the following Hamiltonian with four 
contributions, the Heisenberg and DM interactions, the 
Zeeman term and the uniaxial anisotropy: 


h = - E J iA • $ + E A- • x sj) 

<i,3> <i,3> 

- J2h-Si-Kj2(S?) 2 , 


( 1 ) 


where the interactions are taken with nearest neighbors 
only. 

Since exact diagonalization of this Hamiltonian is not 
possible except for very small systems, we adopt the stan¬ 
dard spin wave approach, that is implemented it two 
clearly separated steps. First, we treat eq. 0 as a clas¬ 
sical functional E c l = H(Qi), where the spin operators S 
are replaced by classical vectors The classical ground 
state, defined as the configuration that minimizes the 
functional £ c i, is determined by self-consistent iteration. 
The above Hamiltonian is studied for square systems with 
up to 45x45 sites. We simulate hard boundary condi¬ 
tions at the edges by the inclusion of a large magnetic 
field in the border sites. Unless otherwise stated, we 
assume [21j [27] S = 1, J = 3.0 meV, D = 0.6 meV, 
Di j || fij and K = 0.5 (D 2 /J), where is the unit 
vector between the sites i and j. We find that, for 
this choice of parameters, the skyrmion configuration is 
more stable than the ferromagnetic state for fields up to 
h c ~ 0 .8(D 2 /J) and skyrmion solutions can not be found 
at all for h c > ~ 1.67(D 2 /J) . Importantly, the main re¬ 
sults discussed below do not depend qualitatively on the 
choice of Hamiltonian parameters. 


(HP) bosons[2l [23]: 

S i • = S — rii, 

S+ = y2 S - Hi Oi, S~ = a\y/2S-ni , (2) 

where Qi is the spin direction of the classical ground state 
on the position i, a\ is a Bosonic creation operator and 
Ui = a\di is the boson number operator. The opera¬ 
tor rii measures the deviation of the system from the 
classical ground state. The essence of the spin wave cal¬ 
culation is to represent the spin operators in Eq. ([l]) 
by the HP bosonic representation and to truncate the 
Hamiltonian up to quadratic order in the bosonic oper¬ 
ators. The linear terms in the bosonic operators vanish 
when the expansion is done around the correct classical 
ground state. This approach has been widely used in the 
calculations of spin waves for ferromagnetic and antifer¬ 
romagnetic ground states [TJ 2 . In the spin wave approx¬ 
imation we thus approximate S^~ = \/2S a\, so that the 
creation of a HP boson is equivalent the removal of one 
unit of spin angular momentum from the classical ground 
state. 

After some algebra [25] the quadratic Hamiltonian can 
be reduced to a symmetric form in creation an annihila¬ 
tion operators written as: 


H — E c i + Eq + 







( 3 ) 


where the constant Eq is a quantum contribution to the 
ground state energy that arises naturally in this expres¬ 
sion as a consequence of applying the bosonic conmuta- 
tion relations in order to bring the Hamiltonian to its 
symmetric form. The second and third term in eq. 0 
correspond to the called spin waves Hamiltonian Hsw • 
Importantly, E c i, Eo, Uj and are a functional of 
the classical ground state configuration that satisfies 

8£cl — Q 
8Qj 

The off-diagonal terms in eq. 0 do not conserve the 
number of HP bosons and are in general present for non- 
collinear classical ground state mm- In consequence, 
the ground state of the Hamiltonian is not the vacuum of 
the a bosons and, thereby, the ground state energy at the 
spin-wave approximation level is different from the classi¬ 
cal ground state energy. The Hamiltonian Hsw is solved 
by means of a paraunitarv[24l [25] transformation, anal¬ 
ogous to the usual Bogoliubov transformation for BCS 
Hamiltonians, that leads to: 


H. 


sw 


N 1 

= Eq T ^ ^ T ~) 

U =1 


( 4 ) 


The second step permits to compute quantum fluctua¬ 
tions of the classical solution. For that matter, we repre¬ 
sent of the spin operators in terms of Holstein Primakoff 


where the a v are bosonic operators formed by linear com¬ 
bination of both a and a t that annihilate spin waves with 
energy huj u . In fig jT] we show the evolution of the four 
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FIG. 1: Left panel: Spin wave spectrum as a function of 
magnetic field. Right panel: Average spin projection along 
the z axis as a function of the external magnetic field. As 
the field grows the size of the skyrmion is reduced. At a 
field of h ~ 1.67 (D 2 / J) the system can no longer sustain the 
skyrmion texture and a transition towards a ferromagnetic 
state is induced. 

lowest energy modes as function of magnetic field , the 
abrupt change at h c > is a consequence of the fact that for 
larger fields, the skyrmion like solution is no longer al¬ 
lowed and a transition to ferromagnetic order takes place. 
It is apparent that the skyrmion has in-gap spin wave ex¬ 
citations, in line with those obtained using the semiclassi- 
cal approximation [21]. The first excitation has a very low 
energy and can be readily associated with a gyrotropic 
traslation mode. In Fig. © we show the magnonic oc¬ 
cupation, defined as 

< = {4> u \a\ai\4> u ) (5) 

for the four lowest energy modes, where \f> v ) = al\GS). 

Inspection of the magnonic occupation in the four low¬ 
est energy excited states, shown in Fig[2] suggest that the 
second and third correspond to a breathing modes. The 
skyrmion-ferromagnet transition is driven by a reduction 
of the energy of the lowest energy breathing mode. The 
abrupt reduction of the second excitation energy near 
the critical field signals the instability of the Skyrmion 
texture. Below the critical field, both internal and ex¬ 
tended spin wave states are present in the system. For 
fields larger than the critical field the spin wave spectrum 
correspond to the usual ferromagnetic spin waves. 

We now discuss the main result of the manuscript. In 
contrast with the semiclassical theory, the quantum the¬ 
ory of spin waves also yields a renormalization of the 
energy of the ground state whose wave function is, by 
definition, the vacuum of the a operators, ol v \GS) = 0. 
In the case of non collinear classical order, the ground 
state is different from the vacuum of HP bosons. The 
ground state energy reads: 

1 N 

£gs = £ c i + £o + x ^ = £d + £zp ( 6 ) 

U=1 

which is different from the classical energy due to the 
contribution coming from quantum fluctuations, or zero 
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FIG. 2: Magnonic occupation of the four lowest energy ex¬ 
cited states displayed in (a), (b), (c), and (d) respectively. 
The plots correspond to a field h ~ 0.75(D 2 /J). On each 
figure the dashed white line correspond to the perimeter of 
the skyrmion. The maximum value of the occupation is 
Umax = 7.1xl0 -3 for the first excited state, Umax — 5.6 xlO -3 
for the second excited state, n maa , = 4.8 x 10 -3 for the third 
excited state, and n macc = 2.1 x 10 -3 for the fourth excited 
state. 


point energy, £zp = £gs ~ £ci • I n the case of ferro¬ 
magnetic solutions it turns out that £zp vanishes mu 
because the collinear ferromagnetic ground state is also 
an eigenstate of the exact Hamiltonian. 

In general the vacuum energy £zp is not zero and, as 
it turns out, it contributes significantly to the relative 
stability of different spin configurations. We show this 
for the case of a skyrmion. In Fig. [3] we compare the 
energy difference between the skyrmion and the collinear 
solution, as a function of the applied field ft, calculated 
at the classical level (blue line) and including £zp (red 
line). It is apparent that the critical field ft c above which 
the collinear solution is more stable increased 10% when 
the £zp is included: the skyrmion is more stable due 
to quantum fluctuations. By analogy with the Casimir 
effect, variations of the vacuum energy with the local 
moment orientation can be associated with an effective 
Casimir magnetic field. This Casimir field is the mag¬ 
netic analogue of the Casimir force. In this manner, the 
effective field on a given spin can be defined as: 

hi = - 5 fff = -ff{£ci + {GS\Usw\GS)) (7) 

oili Dili 

The first contribution is the classical field, which vanishes 
identically by construction. Using this and the Heilman- 
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FIG. 3: Differences between skyrmion and ferromagnetic con¬ 
figuration energy. The crossing between the curves and the 
horizontal axis sets the boundary for the region of stability 
for the skyrmion in classical and corrected approaches (blue 
and red lines respectively). Corrected energy difference be¬ 
tween the skyrmion and ferromagnetic textures including the 
zero point fluctuations, enlarging the region of stability for 
skyrmions. 


Feynman theorem [26], the Casimir field hzp read: 


#ZP = -(GS 


STPrw 

Sfti 


GS) 


(8) 


and they can be readily calculated. Our results for the 
Casimir field are displayed in Fig. @> together with the 
projection of the classical magnetization over the z axis 
(top panel), and the expectation value, computed with 
the ground state wave function, of the magnon occupa¬ 
tion. For reference, we sketch the spherical coordinate 
system as shown in the inset of Fig. @(b). At the cen¬ 
ter of the skymion the magnetization points along — £, 
opposite from the magnetization of the surroundings, as 
can be seen in Fig. 0(a). 

The stabilization of the skyrmion configuration with 
respect to the applied field can be associated to the 
Casimir field using the following argument. The Zeeman 
contribution enters the energy of the system as —Q z h, in 
the same manner, the contribution of the Casimir field in 
classical direction correspond to —S (hzp-&) = — S h^. 
In Fig. 0(c) we plot /iq, it is clear from here that the 
contribution of the Casimir field on the classical direction 
is to lower the energy of the skyrmion state. 

In addition to the stabilization of the skyrmion so¬ 
lution, quantum fluctuations can have other observable 
consequences that we discuss now. In order to charac¬ 
terize the density of spin fluctuations of a given eigen¬ 
state of Hamiltonian 0, |= al\GS), we have de¬ 
fined the magnonic occupation in Ec. 0. Importantly, 
this quantity is non zero even in the ground state, Uq = 


(GS\a\ai\GS}, reflecting the zero point quantum fluctu¬ 
ations that are a consequence of a non collinear classi¬ 
cal state. The zero-point magnonic occupation of the 
skyrmion ground state is shown in figure 0(d). Calcu¬ 
lations using different values for D, J, K show that the 
magnitude of the ZP fluctuations increases as the size 
of the skyrmion is reduced. The ZP fluctuations have 
observable consequences. Since (S • fti) = S — (aja^), 
both the zero point motion and the thermally excited 
spin waves will reduce the effective atomic magnetic mo¬ 
ments of skyrmion solutions, an effect that should be 
observable with local probes, such as STM. 



FIG. 4: Magnetization profile on z direction (a), magnonic oc¬ 
cupation (b) and Casimir field (Q direction) (c) on a skyrmion 
ground state for h=0.75(D 2 /J). The local coordinate system 
for spins and the corresponding Casimir fields is shown in (b) 
inset. As consequence of the Casimir field on direction £4, the 
skyrmion zero point energy is reduced thereby it increase the 
texture stability. 

Importantly, the Casimir magnetic field exert torque 
on the magnetic moments, the Casimir spin torque, de¬ 
fined as Ti = Qi x h l ZP . The ZP energy, and thereby the 
ZP Casimir spin torques, depend on the energy of the 
spin wave excitations, which are interesting on their own 
right and could be observed by means of inelastic electron 
tunneling spectroscopy, as recently demonstrated for spin 
waves in short ferromagnetic chains [28 . 

In summary, we have studied the observable con¬ 
sequences of quantum spin fluctuations in skyrmions. 
We have found that the zero-point energy stabilizes 
skyrmions with respect to the application of a magnetic 
field, increasing significantly the critical field necessary to 
ears the skyrmion. We show that this stabilization can 
be associated to a Casimir magnetic field driven by zero 
point fluctuations analogous to the Casimir effect. We 
have also found that zero point fluctuations renormalize 
the shape of the skyrmion and shrink the atomic mag¬ 
netization, an effect that becomes significant for small 
skyrmions and could be observerved by means of spin- 
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polarized scanning tunneling microscopy. 
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